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Visibility  Probabilities  On  Line  Segments 
In  Three  Dimensional  Spaces  Subject 
To  Random  Poisson  Fields  Of  Obscuring  Spheres.* 

by 

M.  Yadin  and  S.  Zacks 
Technion,  Israel  Institute  of  Technology 
and  State  University  of  New  York  at  Binghamton 


1.  Introduction 


Consider  a  source  of  light  located  at  a  point  O  and  an  arbitrary  straight  line  segment, 
C,  in  a  three  dimensional  Euclidean  space,  .  A  layer  of  random  spheres  is  located 
between  O  and  C .  The  centers  of  the  spheres  are  randomly  dispersed  between  two  parallel 
planes,  according  to  a  Poisson  process  with  a  specified  intensity.  Furthermore,  the  radii 
of  the  spheres  are  independent  identically  distributed  random  variable,  having  a  specified 
distribution,  such  that  both  O ,  and  the  line  segment,  C ,  are  not  covered  by  the  random 
spheres.  A  point  P  on  C  is  called  visible  (in  the  light),  if  the  line  segment  OP  does  not 
intersect  any  random  sphere.  The  measure  of  visibility  on  C  is  the  total  portion  of  C 
which  is  visible.  We  are  interested  in  the  distribution  of  this  random  variable. 

In  a  previous  study  [2,  3,  4]  we  developed  the  methodology  for  determining  the  mo¬ 
ments  of  the  measure  of  visibility  on  star-shaped  curves  in  the  plane  and  provided  also 
an  approximation  to  its  distribution.  The  objective  of  the  present  paper  is  to  apply  the 
two-dimensional  methodology  for  solving  the  three  dimensional  problem  mentioned  above. 
The  results  of  the  present  study  are  thus  restricted  to  the  case  where  the  target  curve 
C,  in  the  three  dimensional  space,  is  a  straight  line  rather  than  an  arbitrary  curve.  A 
runway  in  an  airfield  as  C ,  and  the  pilot’s  position  in  a  flying  airplane  bls  O  ,  is  a  relevant 
example.  In  the  case  under  consideration,  one  can  consider  the  shadowing  process  on  the 
plane  containing  C  and  O.  In  a  previous  Technical  Report  [5]  we  presented  the  basic 
ideas,  but  evaluated  certain  key  functions  by  two  dimensional  numerical  integration.  This 
type  of  numerical  solution  is  time  consuming  and  contributing  to  numerical  error.  In  the 
present  paper  we  provide  an  analytical  solution  to  the  problem  of  determining  simultane¬ 
ous  visibility  probabilities  of  arbitrary  points  on  C.  With  this  anaytical  solution  one  can 
compute  moments  of  random  visibility  measures  on  C  very  fast  and  accurately.  Moreover, 
discrete  approximations  to  the  distributions  of  lengths  of  shadows  on  C  (see  Yadin  and 
Zacks  [6])  can  now  be  computed  efficiently. 


In  Section  2  we  present  the  three  dimensional  shadowing  model,  and  the  basic  geometry. 
In  Section  3  we  discuss  the  reduction  to  a  two-dimensional  shadowing  process.  In  Section 
4  we  present  the  formulae  for  computing  simultaneous  visibility  probabilities  on  a  segment 
of  C.  Analytical  development  of  some  of  the  functions  specified  in  Section  4  is  given  in 
the  Appendix.  In  Section  5  we  provide  an  explicit  solution  for  a  special  case,  in  which  the 
diameters  of  shadowing  spheres,  in  the  Poisson  field  of  obscuring  elements,  are  independent 
and  having  an  identical  tiniform  distribution  on  [0, 6] . 

The  results  of  the  present  study  can  be  applied  to  a  large  class  of  military  and  civil 
problems  of  visibility  in  a  three  dimensional  scenario,  with  randomly  scattered  obscuring 
elements. 

2.  The  Three  Dimensional  Shadowing  Model. 

Consider  a  source  of  light  located  at  a  point  O  (the  origin],  and  an  arbitrary  straight  line, 
C ,  in  a  three-dimensional  Euclidean  space.  A  layer  of  random  spheres  is  located  between 
two  planes  U* ,  IV* ,  which  are  parallel  to  the  line  C.  Let  C*  be  a  plane  containing  C, 
parallel  to  U*  (see  Figure  1).  The  random  spheres,  whose  centers  are  located  between  U* 
and  IV*  cast  shadows  on  C.  As  stated  in  the  Introduction,  the  objective  is  to  develop 
formulae  for  the  computation  of  the  probabilities  of  simultaneous  visibility  of  any  n  points 
on  C  (n  >  1).  These  probabilities  are  then  applied  to  determine  the  moments  of  total 
visibility  measures  on  C ,  and  to  approximate  its  distribution.  In  the  present  section  we 
discuss  the  basic  stochastic  and  geometric  structure  of  the  field. 

Let  u* ,  w*  and  r*  be  the  distances  of  U* ,  W*  and  C*  from  O,  respectively.  It  is 
assumed  that 

(2.1)  0  <  6  <  u*  <  w*  <  r*  -  6, 

where  b  is  the  maximal  radius  of  a  sphere.  Let  M  be  a  plane  passing  through  O  and  C. 
Let  It  and  W  be  straight  lines  at  which  At  intcresects  U*  and  W* .  The  three  parallel 
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lines,  U,  W  and  C  are  at  distances  u,  w  and  r  from  O.  The  following  relationship  is 
satisfied: 


/  X  u*  w*  r*  , 

(2.2)  —  =  —  =  —  =  cos(^), 

uwr 

where  <{>  is  the  angle  between  the  axes  Z*  and  Z,  which  pass  through  O  and  are  per¬ 
pendicular  to  C*  and  C,  respectively.  The  spheres  which  are  scattered  between  U*  and 
W*  have  random  locations  and  random  radii.  It  is  assumed  that  the  centers  of  spheres 
follow  a  random  Poisson  field  of  intensity  A ,  and  that  the  radii  of  disks  are  independent 
random  variables  Xj,  Xj, . . .  having  a  conunong  c.d.f.  F{x) ,  concentrated  on  [0, 6] .  More 
specifically,  let  (ci,C2,C3,x)  be  the  three  location  coordinates  of  the  center  of  a  sphere 
and  its  radius.  Let  S  be  the  sample  space  of  randomly  scattered  spheres,  and  let  B  de¬ 
note  the  Borel  o’ -field  of  subsets  of  $ .  For  J5  G  8 ,  let  N{B}  be  the  number  of  random 
spheres  having  coordinates  in  B.  It  is  assumed  that  for  each  m  >  2,  and  every  parti¬ 
tion  of  S,  N{Bi},...,N{Bm}  are  independent  random  variables,  having 

Poisson  distributions  with  means 

(2.3)  A{Bi)  =  xj  J  dF{x)dci...dc3,  i=l,...,m. 

Bi 

The  random  spheres  in  the  Poisson  field  cast  shadows  on  C*,  and  in  particular  on  the 
line  C .  In  order  to  study  the  shadowing  process  on  C  we  consider  in  the  next  section  the 
projection  of  the  Poisson  field  on  the  plane  Af . 

3.  The  Two-Dimensional  Shadowing  Model  on  M. 

In  the  present  section  we  consider  the  reduction  of  the  shadowing  model  to  Af .  In 
particular,  we  have  to  determine  the  distribution  of  the  centers  of  shadowing  disks  on  M. 
and  the  conditional  distribution  of  their  radii  given  the  location  of  their  centers.  In  Figure 
2  we  present  the  geometry  of  the  reduction  to  a  two-dimensional  model  by  the  projection  of 
the  three-dimensional  space  on  a  plane  perpendicular  to  C .  We  restrict  attention  to  spheres 
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which  are  centered  within  a  prism,  P,  whose  projection  is  the  parallelogram  ABCD  (see 
Figure  2) .  The  distances  from  M  of  spheres  centered  outside  this  prism  is  greater  than  b . 
Accordingly,  such  spheres  do  not  intersect  M  and  cannot  cast  shadows  on  C.  A  sphere 
centered  within  the  prism  P  intersects  M  if  its  radius,  x ,  is  greater  than  its  distance  d 
from  M.  The  intersection  of  such  a  sphere  with  X  is  a  disk  of  radius  y  =  (z^  - 
where  a+  =  max(a,  0) .  For  simplicity,  we  will  assume  that  all  the  spheres  centered  within 
the  prism  P  generate  on  M  disks;  those  spheres  with  x  <  d  generate  on  X  a  disk  with 
radius  y  0.  The  center  of  a  disk  generated  by  a  sphere  is  the  orthogonal  projection  of  its 
center  on  M .  Such  projections  may  lie  between  two  parallel  lines  on  X ,  having  distances 
u  —  /?  and  tw  +  /?  from  O ,  where 

(3-1)  /?  =  6  tan  d>. 


The  orthogonal  projection  of  centers  of  random  spheres  generate  on  X  a  two-dimensional 
Poisson  process.  The  intensity  of  this  process  at  any  given  point  is  proportional  to  the 
length  of  the  segments  within  the  prism  P  on  the  normal  to  X  through  the  point.  As 
seen  in  Figure  2,  this  intensity  decreases  to  zero  towards  the  edges.  In  the  special  case  of 
/?  =  0  the  Poisson  process  on  X  is  homogeneous  with  intensity  rf  =  2A6.  Generally,  the 
intensity  of  the  Poisson  process  on  X ,  at  all  points  on  a  line  parallel  to  U  at  distance  s 
from  O  is 


(3-2)  T}(z)  =  A— [min(2  +  /3,w)  —  max(s  —  /?,u)]. 

Notice  that  ifu-t-^<s<u;  —  ^  then  tj(z)  =  2A6.  That  is,  the  Poisson  field  on  X  is 
homogeneous  only  between  the  parallel  lines  of  distance  u  -f  /?  and  w  - /3  from  O. 

Consider  a  sphere  centered  in  the  prism  P ,  on  a  plane  "V  * ,  whose  distance  from  X  is 
a.  This  sphere  generates  on  X  a  disk  whose  center  falls  on  a  line  of  distance  z  from  O 


(see  Figure  2).  Moreover, 


\z  —  v\  =  aPfb. 


Accordingly,  if  x  is  the  radius  of  the  sphere,  the  radius  of  the  disk  is 


(3.4)  y  =  ((x^-(2-»)"6V/3")+)‘''’- 

Let  G{y  |  z,v)  denote  the  conditional  CDF  of  the  random  radius  K,  of  a  disk  on  A(, 
given  that  its  center  is  on  a  line  of  distance  z  from  O ,  and  the  center  of  the  generating 
disk  is  on  a  plane  intersecting  >1  at  distance  i;  from  O.  According  to  (3.4), 

'0,  y  <  0 

(3.5)  G(y|x,y)=  W  +  (^-'')"4V«'/").  0  <  y  <  6(1  -  ^)^)'/2 

1.  6(1-^)^)'/^  <6. 

It  is  interesting  to  notice  that  G(C  \  z,v)  =  ~  <^|)- 

In  the  following  section  we  show  how  the  simultaneous  visibility  probability  p{si, . . . ,  5„) 
can  be  computed  from  G{y  \  z,v),  and  the  uniform  scattering  in  the  three  dimensional 
prism  P .  For  the  direct  computations  we  partition  the  layer  between  the  two  parallel 
planes  U*  and  If  *  to  parallel  layers  of  width  A*  {  A*  — O).  Consider  a  layer  between  the 
parallel  planes  V*  and  "V**,  within  the  original  layer,  of  distances  v*  and  v*  +  A*  from 
O;  u*  <  V*  <  V*  -b  A*  <  «;* .  Let  v  +  Av  be  the  distance  from  O  at  which  v**  intersects 
M .  We  have 

(3.8)  —  =  — =  —  =  cos  4>. 

V  V  +  Av  r 

Every  point  in  At,  which  lies  on  a  line  parallel  to  V,  whose  distance  from  O  is  between 
V  —  P  +  Av  io  V  +  P  has  the  same  intensity  fi{z)Av,  where 

(3.9)  n{z)  =  Xb/P, 


as  can  be  readily  verified  from  Figure  2.  Since  Av  — >  0,  we  can  consider  the  Poisson  fields 
on  these  thin  layers  to  be  almost  homogeneous. 


4.  Simultaneous  Visibility  of  Points  on  C . 


Consider  n  points  on  C ,  between  P gi  and  P s"  •  Let  si  <  . . .  <  5^  be  the  coordinates  of 
these  points,  with  respect  to  the  orthogonal  projection  of  O  on  C.  In  the  present  section 
we  derive  a  formulae  for  the  simultaneous  visibility  probability,  p(5i, . . ,  ,s„) ,  of  these  n 
points. 

We  say  that  a  specified  set  of  n  points  on  C  are  simultaneously  visible,  if  the  random 
disks  on  M  do  not  intersect  any  one  of  the  n  rays  from  O  to  P,..  (t  =  1, . . ,  ,n) . 

We  introduce  the  points  and  Pg» ,  with  6'  <  s'  and  6"  >  s"  so  that,  disks  centered 
on  the  left  of  the  ray  from  O  to  P$i ,  or  on  the  right  of  the  ray  from  O  to  Pe" ,  will  not 
cast  shadow  on  the  interval  C  =  The  coordinates  of  Pgi  and  Pgn  are: 

(4.1)  e'^s'-bi{s'f  +  ry/^/{u-P), 

and 


(4.2) 


Let  C  be  the  set  of  points  in  A(  such  that,  random  disks  which  are  centered  at  C  can 
cast  shadows  on  the  interval  C,  C  is  a  trapezoid  bounded  by  the  lines  parallel  to  C,  at 
distances  u  ~  P  and  w  +  P  from  O ,  and  by  the  line  segments  OP  ,  and  OP  „  . 

Let  ti{C}  denote  the  expected  numbers  of  disks  having  centers  in  C .  We  consider  also 
disks  of  radius  y  =  0.  According  to  (3.9) 


(4.3) 


w  v  +  /? 

xh  0"-e'  f  f 

• — - - J  dv  I  zdz 


v-p 


= -{0"  -  -  u-^). 

For  values  of  s  in  [O' ,6")  and  0<t<  we  denote  by  XK+{s,t)  [resp.  A/<'_(s,<)) 

the  expected  number  of  disks  centered  between  the  line  segments  OP  and  OP  [resp. 


OP  which  do  not  intersect  OP  . 

- »  —  t  "" a 


O 


Consider  a  disk  centered  on  a  line  "V ,  parallel  to  U ,  which  intersects  the  axis  Z  at 
distance  z  from  O,  u  —  0  <  z  <  w  +  0,  and  having  a  distance  x  from  the  intersection 
of  V  from  OP  (see  Figme  3).  The  maximal  radius  that  such  a  disk  can  have,  without 
intersecting  OP ^ ,  is 

(4.4)  y{x)  =  xr/(5^  +  0  <  x  < 

I 


Hence,  the  fimction  /f+(5,t)  is  given  by, 

to  v+fi  ztjr 

(4.5)  K+{s,t)  =  ^Jdv  j  dz  j  G  I  v,z^  dx. 


u  t)-/9  0 

Geometrical  considerations  yield  that  /^^.(s,^)  =  /f_(5,t)  =  K(x,t) ,  for  every  a  £  [0',0") 
and  all  t  >  0. 

In  the  Appendix  we  prove  that 
b0^t 


(4.6)  K{s,t)  = 


c^r 


^(y,2  _  _  „)^(o)  +  Q  “  Q  ’ 


where  c  =  alJtjb,  a  =  1/(5*  + 


(4.7) 


1 

^'(r)  =  /[I  -  F{bp)]pdp,  0  <  r  <  1, 


and  Q{y,c)  is  given  by  formula  (A. 13).  In  Section  6  we  present  an  explicit  formula  for 
ff(5,  t) ,  for  the  case  of  a  uniform  distribution  of  radii  of  spheres. 

The  simultaneous  visibility  probability  p(si, . . .  ,Sn)  can  be  expressed  in  terms  of  the 
K{s,i)  function,  as  in  Yadin  and  Zacks  [3,  4].  For  n  =  1,  let  p{s)  denote  the  visibility 
probability  of  P, .  This  is  given  by 


(4.8)  p(s)  =  exp{ -[p{C}  -  A(/f(5,5  -  0')  +  K{s,0"  -  a))]  }, 


for  every  0'  <  s  <  0" .  For  n  >  2, 


(4.9) 


p(5i,...,5„)  =  exp{-|r<{C}  -  A(if(5i,si  -  «') 

n— 1 

+  K{sn, 0"  —  a^))  ~  A  ^^(/f(s,-, a,'  —  a,’)  +  /^^(a,’.(.i, a,’.|.i  —  a,'))]  }, 

t=i 
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where 

(4.10)  Si  =  r  tan  +  tan“^  (^*^ '))^  ‘ 

Furthermore,  the  probability  that  the  whole  interval  C  =  \Pt‘,Pa"]  *s  visible  is 

(4.11)  Pi  =  exp{  -X[n*{C}  -  K{s\s'  -  O')  -  A'(s",(9"  -  s")I  }, 


where  is  the  expected  number  of  disks  centered  at  C ,  minus  the  expected  number 


of  disks  centered  between  OP,»  and  OP^n  having  zero  radius,  i.e., 

w  v+P 


=  It{c) j  dv  j 


(4.12) 


U  V  —  P 

w  1 


= /i{C}  —  —  (s"  —  s')  J  dv  j  {u  +  v0)F{b\u\)du. 


-1 


5.  A  Special  Case  And  Numerical  Illustration 

In  the  present  section  we  consider,  for  example,  the  case  of  a  uniform  distribution  of 
sphere  radii,  i.e.,  ^(z)  =  min  ^^,l^/(o,oo(®)»  where  /a(®)  is  the  indicator  function.  In 
the  present  case 

(5.1) 

For  0  <  z  <  X  <  1  define 


1 

*M  =  6“T  +  T’ 


(5.2)  H[..z)  =  j -0=  f  nr)ir. 

X  z/rt 

In  the  present  case  (see  Dwight  [l]) 


(5.3) 


=  ^ I \/l  —  x^  —  ZTT  +  2x  sin  *(x)  +  2z^\/l  —  i^/x 
\/l  -  x2  1  1  +  v^l  -  x2\  I 

i - )]■ 
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Thus,  according  to  (A. 13)  and  (5.2)  the  function  Q[v,c)  is  given  in  the  present  case  by: 


0, 

V  >  yo 

2yo^f  ( 

^yo’ 

yoj 

1. 

yi  <  w  <  yo 

yo  h{ 

<v 

<yo’ 

-1 

yoJ 

V  yi  yo/ 

1  ,  0  <  V  <  yi. 

CtJL 

Substituting  v  =  — ,  j/o  =  (1  +  c^)^/^ ,  j/i/yo  =  c/{l  +  ,  we  obtain  for  a  =  u,\ 


if  o  >  /3 


(l  +  c2)l/2 


2{H-c2)1/2^^. 


.  V^(i  +  c2)^/^’y9(i  +  c2)V2 

{1  +  c^V/Ah( - — _ — _ 

^  ^  L  V/J(1  +  ’  /9(1  +  c2)  1/2 


if/?<a<  ^(l  +  c2)i/2 


^/3(l  +  c2)l/2’^(l  +  c2)V2 

(1  +  02)1/2)]’ 


if  a  <  /?. 


The  function  if(s,  t)  is  given  by  substituting  (5.5)  in  (A.15). 

The  proportional  total  visibility  measure,  on  C  is  defined  as  the  random  integral 


S 

=  jj  Hs)ds, 


where  L  =  s"  —  s'  is  the  length  of  C  and  I[s)  is  a  random  indicator  function,  assuming 
the  value  1  if  and  only  if  P,  is  visible. 

The  moments  of  ly  on  C  are  computed  recursively  according  to  the  algorithm  described 
in  Yadin  and  Zacks  [3].  In  Table  5.1  we  provide  the  first  ten  moments  of  W ,  for  the 
parameters  r*  =  1,  w*  =  .6,  u*  =  .4,  6  =  .3,  s'  =  —1,  s"  =  1.  We  consider  Poisson 
random  fields  with  intensity  parameters  A  =  2,6, 10,  and  inclination  angles  [radians]  4>  = 
0,  f  1  f »  X  ’  nioments  of  the  mixed-beta  approximation  (see  [3])  are  given  also  in  Table 
5.1.  It  should  be  noted  that  the  computation  of  Table  5.1  on  an  IBM-4381  requires  only 
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25  sec'  of  CPU.  An  evaluation  of  the  functions  K{s,t)  by  numerical  integration  raises  the 
computation  time  to  over  3  minutes  of  CPU  and  results  in  loss  of  accuracy. 

In  Table  5.2  we  present  the  parameters  of  the  beta-mixture  approximation  to  the  CDF 
of  W,  namely 


(5.7) 


/ 


0, 


w  <0 


Fw{y>) 


'  Po+  ^  pf°  f  0<u;>l 

B{a,p)  J 


1, 


1  <  u;. 


Using  the  approximation  given  by  (w)  one  can  compute  probabilities  of  W  sets. 


Table  5,1  Moments  of  Visibility  and  Their 
Mixed-Beta  Approximations*,  For  s'  =  — 1,  s"  =  1,  ^  =  0, 
f ,  j,  and  A  =  2,6,10,  r  =  1,  tw  =  .6,  u  =  .4,  6  =  .3 


2 

3 

.932 

.914 

.932 

.914 

.927 

.909 

.927 

.909 

.910 

.887 

.910 

.887 

.847 

.809 

.847 

.809 

.809 

.765 

.809 

.765 

.798 

,751 

.798 

.751 

.754 

.699 

.754 

.699 

.609 

.532 

.609 

.532 

.703 

.640 

.703 

.640 

.687 

.622 

.687 

.622 

.626 

.551 

.626 

.551 

.440 

.352 

.440 

.352 

5 

6 

.894 

.889 

.893 

.886 

.888 

.882 

.886 

.879 

,860 

.853 

.859 

.850 

.765 

.752 

.764 

.749 

.713 

.698 

.712 

.696 

.697 

.681 

.696 

.679 

.635 

.616 

.634 

.614 

.448 

.423 

.447 

.422 

9 

10 

oo 

.887 

.887 

.841 

.873 

.870 

.841 

.879 

.882 

.831 

.865 

,861 

.831 

.846 

.848 

.791 

.832 

.828 

.791 

.733 

.732 

.678 

.721 

.715 

.678 

.676 

.674 

.594 

.664 

.654 

.594 

.657 

.655 

.573 

.646 

.639 

.573 

.586 

.582 

.494 

.576 

.567 

.494 

.382 

.375 

.284 

.375 

.365 

.284 

.515 

.510 

4.20 

.506 

.497 

.420 

.492 

.486 

.396 

.483 

.474 

.396 

.407 

.400 

.309 

.399 

.389 

.309 

.200 

.193 

.123 

.195 

.187 

.123 

*The  upper  value  in  each  cell  is  the  moment  of  W .  The  lower  value  is  that  of  the  beta- 
mixture  distibution,  /loo  =  Pi  • 


Table  5.2  The  Parameters  of  The  Mixed-Beta  Distributions, 
For  s'  =  — 1,  s"  =  1,  ^  =  0,  f,  ^  and  A  =  2,6,10, 
r  =  l,  in  =  .6,  u  =  .4,  6=.3 


Appendix:  Analytic  Derivation  of 

We  derive  here  formula  (4.8)  of  K{s,t).  Starting  from  (4.5),  and  expressing  G{y  |  z,v) 
in  terms  of  F(-) ,  we  obtain 


w  v+P  azt 

(A.l)  K(,.  1)  =  ^  /  iv  /  f  (  (vH  (l^) 

in  which  a  =  (s^  +  r^)"V2 

Let  c  =  a/3t/b  then,  after  several  simple  changes  of  variables,  we  can  write 

2  *  V+c* 

t\  -  Al  f  j..  f  dx  f  F(b(z^  + 


(.4.2) 

c^r  J  “ 

_ 

c^r 

where 

(>1.3) 

1 

J(c,/?,u,u/)  =  f  dy  [ 

1  y+ci 


-1  0 

Notice  that  F{b{x^  +  z^y/^)  =  l  for  all  {x,z)  s.t.  x^  +  z^>l.  Write 


{AA) 


■^(c,/?,  u,t/;)  —  Ji{c,0,u,w)  —  J2{c,  P,u,w), 


where 


F*"  1  Vl-i’ 

(/1. 5)  Ji{c,0,u,w)  =  I  dy  j  dx  j  [1  -  F{b{x^  +  z^y/^)]dz 

fu  -1  0 

and 


1  v^l-x» 

(>1.6)  J2{c,0,u,w)  =  j  dy  j  dx  J  \l  ~  F{b{x^  +  z^y/^)]dx, 

*  M{v  +  CZ,y/l-z^) 


13 


Integrated  Set,  Case  I. 


where  M(a,  6)  =  imn(a,  6) . 

Transforming  (x,z)  to  the  polar  coordinates  (p,0)  we  obtain 


IT  1 


(.1.7) 


=  j  dy  hi  p[l-  F{bp)\dp 


rU  0  0 


=  7r-(u;  -u)'i'(O), 


where  'Ji'(t)  is  defined  in  (4.6). 

In  order  to  evaluate  the  function  J2{c, 0,u^w) ,  define  the  function 


oo  1  oo 


(>1.8) 


Then, 


{A.9) 


«(»,.)  =  J  dy  j  dx  I  |1  -  F(6(i2  + 


V  —1  v+ci 


J2(c,/?,u,u;)  =  _  Q^^u;,c^. 


Consider  the  half  circle  Co  =  {  (x,2)  :  -1  <  x  <  1,  0  <  z  <  1,  x^  +  z^  <1}  (see  Figure 
A.l).  Let  £  be  the  line  z  =  y  +  cx,  having  a  slope  c,  c  >  0,  and  intercept  y.  Let  £o  be  a 
line  parallel  to  £ ,  which  is  tangential  to  Co .  Let  yo  be  the  intercept  of  £o .  Similarly,  let 
be  a  line  parallel  to  £ ,  passing  through  the  point  (-1, 0) .  £i  has  an  intercept  yi  =  c. 
Let  P  be  the  point  at  the  intersection  of  Co  and  £o-  The  right  triangle  with  vertices  O, 
P,  (0,yo)  is  congruent  to  the  triangle  (-1,0),  O,  (0,yi).  Hence  yo  =  (1  +  c^)*/*.  It  is 
clear  that  if  v  >  yo  then  Q{v,c)  =  0.  We  distinguish  two  cases. 

Case  I:  yj  <  v  <  yo. 

The  line  segment  OP  intersects  the  line  £  at  P* ,  whose  distance  from  O  is  7  =  y/yo 
(see  Figure  A.l).  Let  A  and  B  be  the  points  at  which  £  intersects  the  half  circle  Co. 
The  triangle  ABO  A  is  equilateral  and 


(A.IO) 


00  ='^  PQd^  —  *('7)- 


Hence,  by  changing  to  polar  coordinates  and  making  a  proper  rotation,  we  obtain 


yo 


COS  *(-^) 
^  VO  ' 


U-ll)  Q(v,c)=2/dy  f 

J  J  \yo  cos  0 J 

«  0 

Case  II:  0  <  t;  <  j/i . 

Let  denote  the  angle  between  OP  and  OB  (see  Figure  A. 2).  It  is  immediately 
obtained  that  =  f  -  tan-^(c)  =  sin“^  ^  (c)  •  Case  H, 

yo  So  01 

(A.12)  Qiv,c)=  [  dy\ [  - ^)d0 

J  U  \yo  cos  6J  J  \yo  cos  0 J 


Making  the  transformation  fj  —  cos(fl) ,  z  =  y/yoy  >  and  changing  the  order  of  integration, 
we  obtain 


r  o> 


(A.13)  g(u,c)  =  { 


v>yo 


2yo 


i 


yi  <v  <yo 


VO 


yo  / 


r  vdrj 


VO 


VO  *7 
1 


j  ^{z)dz->t-  J  '^{z)dz 


VO*l 


\/l  - 

Notice  that  in  the  case  of  0  <  u  <  y^ ,  one  can  write 

I  ^  i 


,  0  <  r  <  yi. 


(A. 14)  Q{v,c)=yo 


VO 


Vol 


n. 

VO 


Jii_ 

Vo*} 


Finally,  according  to  (A.2),  (A.4),  (A.7)  and  (A.9), 

(A.15)  Kis,t)  =  -u'^)-'^[rv-u)^[0)+Q(^u,^  _  g c)  1 . 
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